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ABSTRACT 

We suggest that extragalactic radio sources are intermittent on timescales of ^ 10^-10^ yr Using a 
simple spherical model of a cocoon/shock system, it is found that inactive sources fade rapidly in radio 
' luminosity but the shock in the ambient medium continues to expand supersonically, thereby keeping the 

0^ , whole source structure intact during the inactive phases. The fading of inactive sources, and the effect of 

the intermittency on the expansion velocity, can readily explain the observed over-abundance of small 
^ ' radio sources. In particular, the plateau in the observed distribution of sizes found by O'Dea & Baum 

(1997) can be interpreted as being due to intermittency. The model predicts that very young sources will 
0^ ' be particularly radio luminous, once the effects of absorption have been accounted for Furthermore, it 

predicts the existence of a significant number of faint 'coasting' sources. These might be detectable in 
deep, low-frequency radio maps, or via the X-ray and optical emission line properties of the shock front. 
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1. Introduction 



Recent radio surveys have identified numerous objects which are morphologically similar to FR-II radio galaxies 
but are appreciably smaller Those sources that are less than 500 pc in extent are termed Compact Symmetric 
Objects (CSOs; Wilkinson et al. 1994). This class also contains many Gigahertz Peaked Sources (GPSs), sources 
^ I whose radio spectrum is seen to peak at GHz frequencies (O'Dea, Baum & Stanghellini 1991). The spectral form 

■ of GPSs is thought to be due to either free-free absorption by an inhomogeneous foreground screen, or synchrotron 

self-absorption in the source itself. Slightly larger sources, those in the range 0.5-15 kpc, have been termed Medium 
KJi Symmetric Objects (MSOs) by Fanti et al. (1995). These various classes of small sources are found to constitute 

j_j ■ 10-30 per cent of all sources in a flux limited sample. 



It is tempting to consider an evolutionary picture in which CSOs evolve into full size FR-II radio galaxies, 
passing through the MSO stage. Since we would expect the sources to remain small for a relatively short period of 
time, there must be strong luminosity evolution for us to see so many small sources (Begelman 1996; Readhead et 
al. 1996). Begelman (1996; hereafter B96) showed how this luminosity evolution could be understood in terms of a 
declining source pressure as the source evolves. 

O'Dea & Baum (1997; hereafter OB97) have recently studied a combined sample of objects including CSOs, 
MSOs and classical FR-II radio galaxies. In particular, they examine the distribution of (projected) linear sizes. 
Using the B96 evolution model, and assuming physically realistic interstellar medium (ISM) density profiles, they 
find a clear overabundance of CSOs and MSOs as compared with the classical FR-II radio galaxies. This is seen as 
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a 'plateau' in the size distribution of their sample of sources between ^ 100 pc and ~ lOkpc. They suggest that 
either (1) a large fraction of the small sources are transient or frustrated and never evolve into large sources, or (2) 
the luminosity evolution is much stronger than that predicted by B96, possibly due to a decline in the efficiency of 
conversion of jet kinetic energy into radio power. 

In this letter we suggest that radio sources are intermittent, and that the source statistics examined by OB97 can 
be understood in the context of a simple evolutionary picture if this intermittency is taken into account. In Section 
2, we develop a simple model of radio source evolution including intermittency. Theoretical source statistics are 
calculated in Section 3. Section 4 discusses predictions of this scenario. Our conclusions are summarized in Section 
5. 



2. The coupled cocoon/shocked-shell model 
2.1. The basic model 

Adopting the standard evolutionary picture (Scheuer 1974; Begehnan & Cioffi 1989), we assume that the radio 
jets are enveloped in, and feed, a cocoon of relativistic material which is overpressured with respect to the ambient 
medium (which may be the ISM of the host galaxy, or the intra-cluster medium (ICM) of the host cluster). This 
overpressure drives a strong shock into the ambient medium and forms a shell of shocked ISM/ICM surrounding the 
relativistic cocoon. The expansion velocity of this shell is determined by the ram pressure of the ambient material 
entering the shock. The cocoon material and the shocked ISM/ICM shell are separated by a contact discontinuity. 

While it is clear that the large scale structures of radio sources show elongation along the jet axis, one rarely 
observes large axial ratios. For example, low-frequency (327 MHz) radio maps of the powerful FR-II source 
Cygnus A (Carilli, Perley & Harris 1994) reveal a radio emitting cocoon with an axial ratio of ~ 3, even though its 
jets are observed to be collimated to within a few degrees. A similar situation is found for the smaller CSOs and 
MSOs. Thus, for the purposes of our simple model, we shall assume that the cocoon and bow-shock are spherical. 
We will denote the radius of the cocoon as r^, and the radius of the bow-shock as rg (rg > rc). We also denote by Vc 
and Vs the volume of the cocoon and shocked shell, respectively. 

We make several further assumptions. Firstly, we assume that at any given instant the pressure within the shock 
(of both the cocoon and the shocked ISM/ICM shell) is spatially uniform with value p{t). See Kaiser & Alexander 
(1997) for an explicit justification of this assumption. Secondly, we suppose that only a small fraction of the total 
kinetic luminosity of the source, L-^{t), is radiated. The rest of this energy is assumed to be fed into the cocoon and 
thus drive the expansion of the cocoon/shocked-shell system. Thirdly, the ambient (undisturbed) medium is assumed 
to have a density distribution of the form p(r) = po{r/a)~", where r is the distance from the center of the radio 
source. 

Given these assumptions, the conservation of energy can be applied to the cocoon and shocked shell to give, 

(7c - ir\V,p + 7cPV;) = Lj{t), (1) 

and 

(7s - l)-i {V,p + -f,pV,) = ^47rrs V(^s)r;^ (2) 

where 7c and 7s are the ratio of specific heat capacities in the cocoon and shocked- shell, respectively, and the dot 
denotes differentiation with respect to time. We close the system of equations with the ram-pressure condition, 
p = p{rs)fs'^. This condition will be valid provided that the expansion of the shocked shell remains highly supersonic 
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Fig. 1. — Properties of a periodically intermittent source (see Section 2.2 of main text for source parameters). Panel 

(a) shows the velocity of the shock in the ambient medium (top) and the radio luminosity of the cocoon (bottom). 
Panel (b) shows the evolution of the cocoon radius (bottom curve) and shock radius (top curve). 

with respect to the ambient medium. If the ambient medium is identified with the hot component of the ISM/ICM, 
its sounds speed wiU be Cg ^ 1000 km . Thus, the expansion of the source will be highly supersonic provided 
r's > few X 1000 km . Once the expansion of the source ceases to be supersonic, the cocoon/shocked-shell 
structure will disrupt and dissipate. 

In order to relate this model to radio observations, we need a prescription relating the radio luminosity, Q, to 
the physical parameters of the model. To do this, we assume that the radio emission is dominated by the synchrotron 
radiation of the relativistic electrons in the cocoon. If we further suppose that the magnetic field is in equipartition 
with the relativistic electrons, standard minimum pressure arguments give Q oc p'^^^Vc. 



2.2. Evolution of a periodically intermittent source 

We now apply the above model to the case of a periodically bursting source. For concreteness, we shall consider 
an initially small source which undergoes 10,000yr long bursts that recur every 100,000yr. We take the jet power to 
be Lj = 10^® crgs^^ during the bursts and negligibly small at other times. The cocoon is assumed to be dominated 
by relativistic material (i.e., 7c = 4/3) and we take a non-relativistic equation of state for the shocked ISM/lCM shell 
(i.e., 7s = 5/3). Finally, the following physically reasonable parameters are taken to characterize the density profile 
of the ambient medium: po = 1-7 x 10~^^ gcm~^ , a = 500 pc and a = 1.5. The quaUtative behaviour described 
below is insensitive to reasonable departures from these canonical parameter values. 

Figure 1 shows the results of a numerical integration of eqns. (l)-(2). During the initial burst of activity 
(t < 10, 000 yr) the source expands self-similarly, i.e., the cocoon radius is a constant multiple of the shock radius. 
This is the evolutionary phase that has been previously examined by Falle (1991), B96, and Kaiser & Alexander 
(1997). As found in previous works, Vs oc By assumption, the power source switches off at t = 10, 000 yr 

and the cocoon/shocked-shell system enters a 'coasting' phase akin to Taylor-Sedov expansion. It must be stressed 
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that the expansion during this coasting phase is still pressure-driven, as opposed to being driven by the momentum of 
the shocked shell. Despite the increased deceleration of the shocked shell, it still remains highly supersonic for the 
entire period of this coasting phase (i.e., until t = 100, 000 yr). However, due to the drop in source pressure, the radio 
luminosity of the cocoon will fall rapidly once the power source has turned off. Thus, even though the basic source 
structure remains intact, a coasting source will be significantly fainter than an active source. It is also interesting to 
note that the source evolution is no longer self-similar in the coasting phase. 

The onset of the second burst of activity induces a rapid increase in pressure leading to a subsequent increase 
in both the source expansion rate and cocoon radio luminosity. Note that during the (brief) period in which the 
contact discontinuity is accelerating (i.e., f'c > 0), this interface will be subject to Rayleigh-Taylor instabiUties that 
will tend to mix material from the shocked shell with the relativistic cocoon material. This mixing will produce 
time-dependence in the effective value of 7c, thereby affecting the detailed evolution of the cocoon. The implications 
of this mixing will be addressed in future work. 

After many bursts (or, more precisely, when the recurrence timescale is short compared with the expansion 
timescale), the intermittent nature of the source will be unimportant in determining its evolution. It will behave as a 
constantly fed source with jet power /ij, where / is the fraction of time that the source is on with power Lj. The 
re-estabUshment of an approximately self-similar expansion in our numerical experiment reflects this fact. 



3. Source statistics 

Suppose that we have a single population of evolving radio sources. Further, suppose we form a flux-Umited 

sample of these sources. Assuming a Euclidean universe, the number of sources in the luminosity range Q ^ Q-\- dQ 
is dN oc <t){Q)Q^/'^dQ, where <j){Q)dQ is the volume density of sources in that luminosity range. Since (j){Q)dQ is 
proportional to the time that a given source spends in this luminosity range, we have 

m) oc = (^r-e^) , (3) 

implying that 

^ocQ^/Vc-^ (4) 
arc 

The evolutionary model of Section 2 allows us to determine the functions Q{rc) and rc(rc). Thus, eqn. (4) is an 
explicit expression for dN/dr^ as a function of r^. We have chosen to examine the distribution of since it is the 
cocoon radius that will be identified observationally as the half-size of the radio source. 

Figure 2a shows a comparison of the observed size distribution of OB97 with our theoretical size distribution for 
a single population of periodically intermittent sources. To facilitate this comparison, we have biimed the theoretical 

size distribution using bins of 0.5 log(2rc). In order to match the distribution of OB97, we set a = 1.8 (determined 
by the slope of the distribution at large sizes) and assume a burst duration of 30,000 yr. All other parameters have the 
values of Section 2.2. 

The intermittency of the sources allows the qualitative features of the OB97 size distribution to be reproduced. In 
particular, there is a plateau in the size distribution resulting from sources that are still undergoing their first few bursts 
of activity. If the break in the OB97 distribution at small sizes (< 100 pc) is real, this could be identified as being due 
to sources that are still undergoing their first burst of activity. In this ideahzed case of a single evolving population, 
there is fine structure within the size distribution corresponding to the distinct cycles of activity. In practice, the 
stochastic nature of the parameters in any real source population will wash out this fine structure. In particular, the 
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Fig. 2. — (a) Theoretical and observed size distributions. The dotted line shows the fine-grained (i.e. unbinned) 
theoretical model which, for clarity, has not been displayed beyond lOkpc. The filled squares show the binned 
theoretical distribution using bins of size 0.5 log(2rc). The data from Fig. 10 of OB97 (with 1-cr errors) are shown as 
open squares, (b) radio luminosity as a function of rc for these same (theoretical) intermittent sources. 

disagreement between our model and the data at ~ 2 kpc can be resolved if we consider realistic source populations. 
Note that we have not included the largest size bin of the OB97 distribution since this is probably affected by the 
complete turning-off of old sources. 

Figure 2b shows the radio luminosity Q as a function of the total source size 2rc- This is to be compared with 
Fig. 9 of OB97. An important feature of Fig. 2b is the dramatic decline in radio luminosity between the first burst 
and all subsequent bursts. In other words, if one were to assume a constantly fed source and extrapolate from large 
sources to small sources, then one would substantially underestimate the small source luminosity. 



4. Discussion 

There are two immediate predictions of the intermittency model. First, Fig. 2b suggests that statistically- 
complete samples of GPSs and small CSOs will show them to be very much more radio luminous than suggested 
by extrapolation of MSO/FR-11 properties. Assuming that the (optically thin) radio luminosity depends on cocoon 
pressure in some way, this prediction should hold for any physically reasonable prescription relating p, Vc and the 
radio luminosity, Q. It must be noted, however, that there are observational complications involved in testing this 
prediction. As well as being comparatively rare, these small sources are often absorbed at typical radio frequencies. 
This may be due to free-free absorption in an inhomogeneous foreground screen or synchrotron self-absorption in the 
source itself. This absorption must be corrected for before comparisons of the type discussed above can be made. 

A second prediction is that there should be a large number of medium size objects (a few hundred to a few 
thousand parsecs across) which are in a 'coasting' phase and have faded below the flux limits of current radio surveys. 
There are several methods that could be employed to search for such sources. Deep, low-frequency radio maps might 
reveal the coasting cocoons of such sources. Alternatively, we might hope to observe the ISM/ICM shock either 



-6- 



through X-ray signatures (using the high spatial resolution of AXAF) or via the Ha emission that it surely excites 
(e.g., see Bicknell & Begelman 1996). 



5. Conclusions 

We have explored the implications of radio source intermittency on source statistics. To do this, we have 
developed a simple model for the evolution of a cocoon/shocked-sheU system which is expanding supersonically 
into an ambient medium that possesses a power-law density profile. The cocoon is assumed to be fed energy at a 
rate Lj{t). This model is integrated numerically for the case of a periodic source which has active phases (with 
constant Lj) separated by inactive, or coasting, phases in which Lj = 0. During the first few periods of inactivity, the 
radio luminosity will fade rapidly. However, these young sources can maintain highly supersonic expansion during 
their coasting phases and, hence, will remain intact throughout the inactive periods. Once a source has grown large 
enough such that the expansion timescale is longer than the recurrence timescale, the intermittency will not affect its 
subsequent evolution. The fading of small, inactive sources, and the effect of intermittency on the expansion velocity 
of the sources, produce a double-break in the size distribution. This can be readily identified with the plateau found in 
the size distribution of OB97. 

There are two clear predictions of this model. First, one could search for the faint, inactive sources. These 
sources might reveal themselves in deep, low-frequency radio surveys. Alternatively, they could be detectable via the 
X-ray signatures or Ha emission accompanying the coasting ISM/ICM shock front. Secondly, statistically complete 
radio samples of small sources (~ 100 pc and smaller), once corrected for absorption effects, should show these 
sources to be very overluminous as compared with an extrapolation from larger sources. Since this is essentially 
reflecting the high pressure of these small sources, this prediction should be independent of the precise form of the 
radio emissivity. 
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